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We report results of air monitoring started due to the recent natural catastrophe on 11 March 2011 in Japan and the 
severe ensuing damage to the Fukushima Dai-ichi nuclear reactor complex. On 17-18 March 2011, we registered the 
first arrival of the airborne fission products 131 I, 132 I, 132 Te, 134 Cs, and 137 Cs in Seattle, WA, USA, by identifying their 
characteristic gamma rays using a germanium detector. We measured the evolution of the activities over a period of 23 
days at the end of which the activities had mostly fallen below our detection limit. The highest detected activity from 
radionuclides attached to particulate matter amounted to 4.4±1.3 mBqm~ 3 of 131 I on 19-20 March. 

Keywords: radiation monitoring, fission products, germanium detector 
PACS: 89.60.Gg, 25.85.-w, 29.40.-n 
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1. Introduction 

The recent earthquake and tsunami in Japan on 11 
March 2011, 05:46 Coordinated Universal Time (UTC) re- 
sulted in severe damage to the Fukushima Dai-ichi nuclear 
reactor complex. Due to the uncertainty of the situation, 
limited quantitative information, and its potential impact 
on both local public health as well as our low-background 
fundamental physics program (Elliott et al. 2009), we be- 



gan monitoring local air samples in Seattle, WA, USA, for 
the potential arrival of airborne radioactive fission prod- 
ucts. In this paper we present data on key radionuclides 
associated with the nuclear accident and a brief discussion 
of the transport. We believe it is important to provide a 
rapid release of this data for two reasons: 

• While the earthquake, tsunami and nuclear accident 
are extremely serious within Japan, there has also 
been a great deal of concern around nuclear radia- 
tion reaching the U.S. Our data, generated by inde- 
pendent university research, should provide greater 
confidence in our understanding of the risks from 
long-distance transport of radionuclides. 

• These data, along with other observations, will help 
provide better information on the characteristics of 
the source, the release mechanism, and the transport 
of the radionuclides. 
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2. Material and Methods 

Our samples consist of air filters taken from the intake 
to the ventilation system of the Physics and Astronomy 
building at the University of Washington. This allows us 
to sample ~10 times more air than what had been done 
previously here at the university after the Chernobyl inci- 



dent (Kellogg et al. 1986 ) and proved to be one of the key 



points for the successful detection of the radioactive fission 
products. In order to search for characteristic gamma rays 
stemming from radionuclides we place the samples inside 
a lead shield of 5 to 20 cm thickness next to a 0.5-kg p- 



type point contact germanium detector (Barbeau et al. 
|2007[ ) for low-level counting. The detector exhibits an en- 
ergy resolution of 1.4 keV FWHM at 600 keV. The level 
of observed background radiation inside the shield ranges 
from 10 counts/keV/hour at 50 keV to 2 counts/keV/hour 
at 800 keV. The energy and efficiency of the detector have 
been calibrated using 10 strong gamma lines between 200 
and 1500 keV from a 152 Eu source to an accuracy of about 
0.1 keV and 10%, respectively. 

We digitize the preamplified traces coming from the 
germanium detector using the Struck card SIS3302 ( |Struck| 



Innovative Systeme 2011 ) which at the same time extracts 
the energy of the measured pulse. The communication 
with the card and the VME crate is managed by ORCA 



(Howe et al. 2004). After the acquisition the data are au- 



tomatically uploaded to a database for analysis. In addi- 
tion to the real physics events we also inject a pulser signal 
into the preamplifier to check the live time and health of 
the system. The pulser runs at 0.1 Hz with an amplitude 
equivalent to an energy of ~120 keV. The signal from the 
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pulser is clearly visible in Fig. [T] 

The air niters used are commercial ventilation filters 
from AmericanAirFilter (Model PerfectPleat ULTRA) and 
Purolator (Model DMK80-STD2) with dimensions 61x61x5 
Their efficiency for retaining particles down to a size of 
5 /im amounts to 75%, drops to 35% at 1 /im and to 5% at 
0.4 /im. From our detection of the cosmogenic 7 Be isotope, 
we calculated an activity of ~0.1 mBqm~ 3 . Comparing 
this value to the known 7 Be concentration of 2-8 mBqm~ 3 



(Yoshimori et al. 2003) we deduce a filter efficiency of 



2.4±2.0% or, correspondingly, particle sizes of <0.4 /xm. 
This roughly agrees with observations of radioactive parti- 



cle sizes after the Chernobyl accident (Devell et al. 1986 



Jost et al. 1986 ) and measured sizes of radioactive dust 



in the atmosphere ( Bergamini et al. 1973 ) . In an auxil 



iary measurement we exposed in addition to our standard 
filter a high efficiency particulate air filter (Hunter HEP- 
Atech filter model 30930). This filter is -100% efficient 
for particles of sizes >0.1 /xm and features one layer of 
activated charcoal. Taking the measured activities from 
that filter as the true activities in air, we calculate the ef- 
ficiency for our standard filters from the ratio of measured 
activities. The filter efficiency amounts to 14.2±3.2% for 
the 131 I and 4.4±1.2% for the other fission products and 
7 Be, in agreement with the above estimate. This points to 
slightly larger particle sizes of 131 1 compared to the other 
isotopes (~0.6 /im instead of ~0.4 zxm given the filter effi- 
ciencies). We can only speculate on the initial conditions 
and the chemistry involved that would lead to this differ- 
ence. A difference between the particle sizes for 131 I and 
other fission products had already been observed after the 



Chernobyl accident (Jost et al. 1986). However, in those 
measurements the 13i I particle sizes tended to be slightly 
smaller than the other fission products. We believe that 
the difference between those measurements and our find- 
ings is due to the longer transport distance. The long 
transport tends to lead to overall smaller particle sizes. 
As 131 I can be transported part of the way in its elemen- 
tal form and can thus be adsorbed at later stages, this 
could lead to larger 131 I particle sizes with respect to the 
other fission products which are adsorbed shortly after the 
release. 

At several other radiation monitoring stations cartridges 
filled with charcoal were used in addition to particle air fil- 




ters ( RadNet| [201 1| |Taylor| |2011[ ). The charcoal captures 
the 131 1 present in its gaseous form and the combination of 
the two thus measures the total amount of 131 I in air. The 
ratio of gaseous 131 I to 131 I attached to particulate matter 
varied between 2 and 20 with an average of about 5. We 
would like to stress that our results below only measure 
the portion of 131 I attached to particulate matter. 

The air filters were typically exposed for one day to 
an air flow of 114000±8000 m 3 /day, which was measured 
using a Davis 271 Turbo-Meter flowmeter. We bagged and 
compressed the filters into packages of approximately 1000 
to 4000 cm 3 before placing them into the lead shield for 
counting. The solid angle for gamma rays emitted within 
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Figure 1: (Colour online) Comparison of the gamma spectra from 
the measurements of air filter PHI (red, 16-17 March) and air filter 
PH2 (blue, 17-18 March) showing clearly the additional peaks due 
to the arrival of radioactive fission products at the US west coast. 
The dominant peak at 364 keV is from 131 1. 



that volume and interacting with the germanium detec- 
tor was calculated given the actual dimensions and ranged 
between 1.6 and 4.4%. We attribute a 20% systematic 
uncertainty to the calculated value. 

3. Results and Discussion 

We started the air monitoring campaign on 16 March 
2011. The exact exposure and counting periods for the dif- 
ferent air filters are listed in Table [l] No fission products 
were detected in the first air filter PHI and we were able to 
attribute all the visible gamma lines to known background 
radioactivity from cosmic-ray induced processes, various 
radioactive isotopes of the uranium and thorium decay 
chains, cosmogenic 7 Be, and 40 K. The subsequent sample 
PH2 immediately revealed the onset of several character- 
istic gamma lines from fission products. The identified 
isotopes are 131 I, 132 I, 132 Te, 134 Cs, and 137 Cs. Figure [l] 
shows the comparison between the gamma ray spectra 
from air filter PHI and PH2 where the additional gamma 
peaks are clearly identifiable. Figure [2] demonstrates the 
statistical significance of the detected lines by showing the 
peaks of the strongest decay branches of the five identified 
isotopes. 

In order to obtain the decay rates of the different iso- 
topes we summed the appropriate spectral bins in the re- 
gion of interest and subtracted the sum of the same amount 
of bins in the side bands. For the cases in which the sta- 
tistical significance of the extracted signal counts was less 
than three sigma we calculated the upper limit at 95% 



C.L. by employing the Feldman- Cousins formalism (Feld 



man and Cousins 19981 for a possible signal given the 
number of background counts. From the detected number 
of counts (iVdet ) and the given counting time (T) and iso- 
tope lifetime (t) we calculate the decay rate (i?dct) at the 
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Figure 2: (Colour online) Plot of the 5 strongest gamma lines of 
131j i 132 I; 132 T6] 134 CSj an( j 137 Cg for the air filter p R1 ( red ) and 

air filter PH4 (blue) measurements. The change in 214 Biactivity is 
due to fluctuating radon levels during the time of measurement. 



beginning of the counting period: 

N dct 



Rdc 



1 - e~ T /r) 



(1) 



From this and measurements and calculations of air flow 
(Q), filter, geometrical, and detection efficiency (e n iter, 
and edet) the amount of activity present in air (A a \ T ) 



c gcom 



using the known branching ratios (BR) is: 

Rdct 



A, 



1 



L ! < Mil 



(2) 



In addition, we correct the values for delays in the count- 
ing period (idol) with respect to the end of the exposure 
using the known lifetimes (r) thus calculating the activity 
present in the filter at the end of the exposure time: 



4' 



(3) 



We do not correct for decays during the exposure itself as 
the time structure of the arrival of the radioactive atoms 
is unknown. The values for the activities can be found 
in Fig. [3] and Table [2j We cross-checked for consistency 
between the activities obtained from different branches of 
the same isotope. As they were in agreement and the 
systematic error dominates we only give the value for the 
strongest branch. For 132 I with a half life of only 2.3 hours 
we do not detect the atoms originating in the reactor but 
the ones following the decay of 132 Te. The measured ac- 
tivities for these two isotopes given in Table [2] agree very 
nicely corroborating our correction due to the energy de- 
pendent detection efficiency. The highest observed activity 
of 131 I amounts to 4.4±1.3 mBqm" 3 . 

We looked for a correlation between the activities ob- 
tained for the fission products and the cosmogenic isotope 
7 Be in order to probe the influence of weather effects on 



"E 2.0 

CD 
E, 

£ 1.5 

> 
o 

5 1.0 



0.5 



0.0 



I (1/3 x) 

12| 

,2 Te 
,7 Cs 
,4 Cs 



|H^^\>. tty , „ 



10 15 20 25 30 

Time since earthquake [days] 



Figure 3: (Colour online) Activities of the five fission isotopes 131 I, 
132 I, 132 Te, 134 Cs, and 137 Cs in air for the exposed air filters covering 
a period of 23 days. Results for 131 1 are scaled by a factor 1/3 in 
order to improve the visibility of the other radionuclides. The exact 
values can be found in Tabled The errors show the statistical and 
systematic uncertainties added in quadrature. Points are shown with 
small horizontal offsets to increase visibility. 



the observed fluctuations in activities. For this we calcu- 
lated the Pearson product-moment correlation coefficient 

(see, 



e.g. 



Cowan (1998)) between the fission isotope ac- 
Y Be. 



tivities and Be. The values obtained for the correlation 
coefficient between 0.34 and 0.61 indicate a "medium" cor- 
relation. This suggests that part of the observed fluctua- 
tions are due to changing weather conditions that lead to 
more or less activity passing through the air filters. 

In addition to the activities themselves we also ana- 
lyzed the ratio of activities for several isotopes. While the 
activities have a large systematic error due to the different 
efficiency corrections these errors cancel to a very large 
degree in the ratios. Figure H shows the ratios for 131 I / 
137 Cs, 132 Te / 131 I, 132 Te / lW Cs, and 134 Cs / 137 Cs. We 
show also fits to the ratios r with a simple exponential 
decay 



r(t) = r e 



(4) 



where the lifetimes Tj have been fixed to the known values 
of the different isotopes. We excluded the points where 
only an upper limit for one of the activities was available. 
While the ratio of the two caesium isotopes is nicely de- 
scribed by this fit (a constant fit in that case due to the 
large half lives with x 2 /dof = 10.1/8) the other ratios show 
large deviations (x 2 /dof >5) and the fit rather follows the 
general trend of the ratio. This indicates that the release 
and transport efficiencies for the different elements vary 
in time leading to "bursts" of individual elements dur- 
ing certain days. It has already been alluded above that 
certainly the transport of 131 I is different from the other 
radionuclides due to the differences in the measured filter 



3 



Table 1: Exposure and counting periods for the different air filters. Times are given in UTC. Days after earthquake give the time at the end 
of the exposure. 



Filter 


Brand 




Exposure 






Days after earthquake 




Counting 




PHI 


AmericanAir 


16/3 


17:00 


17/3 


19 


00 


6.6 


17/3 


23:45 


- 18/3 


22:45 


PH2 


AmericanAir 


17/3 


19:00 


- 18/3 


21 


00 


7.6 


18/3 


23:30 


- 19/3 


02:30 


















19/3 


06:00 


- 19/3 


16:00 


PH3 


AmericanAir 


18/3 


21:00 


- 19/3 


16 


20 


8.4 


19/3 


18:30 


- 20/3 


18:30 


PH4 


Purolator 


19/3 


16:20 


- 20/3 


18 


00 


9.5 


20/3 


19:30 


- 21/3 


16:30 


PH5 


Purolator 


20/3 


18:00 


- 21/3 


21 


00 


10.6 


21/3 


23:20 


- 22/3 


21:05 


PH6 


Purolator 


21/3 


21:00 


- 22/3 


20 


30 


11.6 


22/3 


21:16 


- 23/3 


02:16 


















24/3 


01:44 


- 24/3 


20:45 


PH7 


Purolator 


22/3 


20:30 


- 24/3 


01 


00 


12.8 


24/3 


22:48 


- 25/3 


21:48 


PH8 


Purolator 


24/3 


01:00 


- 25/3 


00 


15 


13.8 


27/3 


05:08 


- 27/3 


20:08 


PH9 


Purolator 


25/3 


00:15 


- 25/3 


21 


00 


14.6 


26/3 


17:59 


- 27/3 


04:59 


PH10 


Purolator 


25/3 


21:00 


- 27/3 


04 


45 


16.0 


28/3 


16:11 


- 29/3 


01:11 


PH11 


Purolator 


27/3 


04:45 


- 28/3 


15 


05 


17.4 


29/3 


01:48 


- 29/3 


15:48 


PH12 


Purolator 


28/3 


15:05 


- 29/3 


16 


15 


18.4 


29/3 


16:59 


- 30/3 


19:00 


PH13 


Purolator 


29/3 


16:15 


- 30/3 


17 


00 


19.5 


30/3 


20:27 


- 01/4 


02:27 


PH14 


Purolator 


30/3 


17:30 


- 31/3 


20 


00 


20.6 


01/4 


18:14 


- 03/4 


01:14 


PH15 


Purolator 


31/3 


20:00 


- 01/4 


21 


00 


21.6 


04/4 


15:37 


- 05/4 


18:14 


PH16 


Purolator 


01/4 


21:00 


- 04/4 


23 


00 


24.7 


05/4 


19:05 


- 07/4 


17:05 


PH17 


Purolator 


04/4 


23:00 


- 06/4 


21 


30 


26.7 


07/4 


20:22 


- 08/4 


19:22 


PH18 


Purolator 


06/4 


21:30 


- 08/4 


21 


00 


28.6 


10/4 


03:35 


11/4 


22:35 



efficiencies. The release and transport of 131 1 is also more 
efficient than that of 132 Te as can be seen by extrapolating 
the ratio of 132 Te to 131 I back to the end of active reactions 
in the fuel rod (day in Fig.|4|. Comparing the extracted 
value of ~0.5 to the ratio of fission yields of ~ 1.5 (Isotope 



Project, 2011) shows that more I was released and/or 



transported than 132 Te. 

In order to assess the transport time of the fission 
products across the Pacific we performed several model 
calculations. The trajectories were computed using the 



NOAA HYSPLIT model ( praxler and Rolph| [20TT| ) us- 
ing the Global Data Assimilation System (GDAS) meteo- 
rological dataset and model-calculated vertical velocities. 
We used the HYSPLIT model to calculate several hun- 
dred trajectories over the time period of interest. The 
trajectories show a range in transport patterns depending 
on the altitude and hour of the start time. Figure [5] (a) 
shows three trajectories which exhibit the range of trans- 
port pathways. The start time is 12 March 2011, 10 UTC, 
which is approximately 3 hours after the first reported 



explosion from Unit 1 (Tokyo Electric Power Company 



TEPCO 2011), and the trajectories were calculated for 



three heights in the boundary layer, 500, 1000 and 1500 
meters above ground level. The trajectory at 500 meters 
is caught up in, and lofted by, a cyclonic system over the 
Bering Sea. The trajectories started at 1000 and 1500 
meters are partially lofted by the same system, but do 
not enter the cyclonic pattern. Instead they are rapidly 
transported across the Pacific. Upon arrival to the west 
coast of the U.S. the transport again splits with one arm 



transported to the north in a cyclonic direction around a 
second low pressure system located off the coast of Wash- 
ington state. Figure [5] (b) shows the 850 mb geopotential 
heights, wind vectors and temperatures for 17 March 2011, 
12 UTC. The trajectory initially started at 1500 meters is 
transported in the boundary layer towards California. The 
northern arm is again lofted by a cold front near Wash- 
ington state. There were rain showers and cool weather 
in Western Washington at this time. The strong diver- 
gence and precipitation associated with these weather sys- 
tems most likely significantly reduced the concentrations 
of radionuclides that were transported. The trajectories 
support the notion of transport of the radionuclides from 
the Japanese boundary layer to the U.S. boundary layer 
in only 5 to 6 days. This is significantly faster than pre- 



viously reported trans-Pacific times ( Weiss-Penzias et al. 



2007 National Research Council-Board on Atmospheric 



Sciences and Climate 2009 1 , especially considering the ra- 



dionuclides were released in the boundary layer over Japan 
and measured in the boundary layer along the U.S. west 
coast. 

Following the nuclear accident, measurements on radi- 
ation levels on land and off the coast of Fukushima from 
several sources were collected and made available by the 



( Ministr y of Education, Culture, Sports, Science and Tech- 
Inology MEXTl 12011). Highest reported readings for 131 I 



in the region area are in the range of 500-600 Bqm on 
March 22nd and 23rd about 25 km south of the plant, 
although there is an extreme value of 5600 Bqm~ 3 re- 
ported for March 21st. The Tokyo Electric Power Com- 



4 



Table 2: Detected activities and upper limits for the measured radionuclides attached to particulate matter. The numbers are given in /^Bq 
per m 3 . The errors give the 1<t values for the systematic and statistical uncertainties (in that order). In case of upper limits we give the limit 
from statistics only at 95% C.L. and quote the influence of the systematic error separately at la. The upper limits strongly depend on the 
fluctuations of background radiation at the time of counting. 



7 BS ^ re*! ^Ti ™l ^ci n?7 Ci 

Sample 53.2 days 8.0 days 2.3 hours 3.2 days 20.8 hours 2.1 years 30.1 years 

477.6 keV 364.5 keV 667.7 keV 228.2 keV 529.9 keV 604.7 kcV 661.7 keV 



PHI 


1200 


± 


300 ± 


200 


< 2.5 


± 


0.7 


< 16 


± 


5 


< 6 


± 


2 


< 70 


± 


20 


< 17 


± 


5 




< 11 


± 


3 




PH2 


2200 


± 


600 ± 


300 


1690 


± 


490 ± 20 


550 


± 


150 ± 30 


580 


± 


160 ± 50 


< 170 


± 


50 


< 80 


± 


20 




140 


± 


40 


± 40 


PH3 


4100 


± 


1100 ± 300 


2160 


± 


630 ± 20 


640 


± 


180 ± 30 


530 


± 


150 ± 40 


< 110 


± 


30 


160 


± 


40 ± 


30 


210 


± 


60 


± 40 


PH4 


1400 


± 


400 


± 


200 


4440 


± 


1300 ± 30 


790 


± 


220 ± 30 


780 


± 


220 ± 40 


< 60 


± 


20 


110 


± 


30 ± 


30 


140 


± 


40 


± 30 


PH5 


1100 


± 


300 


± 


200 


1320 


± 


380 ± 20 


130 


± 


40 ± 20 


130 


± 


40 ± 30 


< 40 


± 


10 


70 


± 


20 ± 


20 


90 


± 


20 


± 30 


PH6 


1200 


± 


300 


± 


200 


980 


± 


290 ± 10 


60 


± 


20 ± 20 


< 80 


± 


20 


< 9 


± 


3 


80 


± 


20 ± 


30 


140 


± 


40 


± 30 


PH7 


1560 


± 


430 


± 


100 


1190 


± 


350 ± 10 


< 50 


± 


10 


< 70 


± 


20 


< 180 


± 


50 


39 


± 


11 ± 


9 


< 50 


± 


10 




PH8 


1900 


± 


500 


± 


100 


1550 


± 


450 ± 20 


290 


± 


80 ± 20 


220 


± 


60 ± 30 


< 170 


± 


50 


220 


± 


60 ± 


20 


220 


± 


60 


± 20 


PH9 


2100 


± 


600 


± 


200 


1260 


± 


370 ± 20 


70 


± 


20 ± 10 


80 


± 


20 ± 20 


< 80 


± 


20 


110 


± 


30 ± 


20 


150 


± 


40 


± 20 


PH10 


1900 


± 


500 


± 


100 


800 


± 


230 ± 10 


50 


± 


10 ± 10 


< 60 


± 


20 


< 32 


± 


9 


120 


± 


30 ± 


10 


150 


± 


40 


± 20 


PH11 


1070 


± 


290 


± 


90 


431 


± 


126 ± 7 


< 10 


± 


3 


< 9 


± 


3 


< 11 


± 


3 


< 16 


± 


4 




< 20 


± 


5 




PH12 


1070 


± 


300 


± 


100 


470 


± 


139 ± 7 


< 35 


± 


10 


< 13 


± 


4 


< 9 


± 


3 


< 0.8 


± 


0.2 




< 13 


± 


4 




PH13 


1910 


± 


530 


± 


100 


362 


± 


107 ± 6 


< 16 


± 


5 


< 11 


± 


3 


< 6 


± 


2 


< 40 


± 


10 




< 40 


± 


10 




PH14 


410 


± 


110 


± 


60 


26 


± 


8 ± 2 


< 10 


± 


3 


< 19 


± 


5 


< 40 


± 


10 


< 10 


± 


3 




< 23 


± 


6 




PH15 


730 


± 


210 


± 


80 


43 


± 


13 ± 3 


< 10 


± 


3 


< 40 


± 


10 


< 270 


± 


80 


< 12 


± 


3 




< 12 


± 


3 




PH16 


840 


± 


240 


± 


30 


227 


± 


67 ± 2 


< 1.3 


± 


0.4 


< 6 


± 


2 


< 11 


± 


3 


17 


± 


5 ± 3 




31 


± 


9 ± 4 


PH17 


700 


± 


200 


± 


60 


89 


± 


27 ± 3 


< 7 


± 


2 


< 7 


± 


2 


< 22 


± 


6 


23 


± 


7 ± 7 




< 26 


± 


7 




PH18 


460 


± 


130 


± 


40 


44 


± 


13 ± 2 


< 2.4 


± 


0.7 


< 2.9 


± 


0.8 


< 23 


± 


7 


< 19 


± 


5 




18 


± 


5 ± 5 
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Figure 4: (Colour online) Ratio of the activities for the isotopes 
specified in the legend. The error band shows the statistical error 
only as the systematic errors cancel to a large degree in the ratios. 
The fits to the ratios as a function of time were performed with fixed 
decay time constants. For a discussion see text. 



pany (TEPCO) is constantly monitoring 131 I and reported 
one reading of 4100 Bqnr 3 within about 10 km of the re- 
actor on March 20th. Most other observations by MEXT 
and TEPCO are less than 10 Bqm~ 3 . The Japanese Min- 
istry of Defense also reported 131 I values using aircraft 
samples, with a maximum of 0.46 Bqm~ 3 for March 25th 
in an area 25-30 km to the west. Based on our highest ob- 
served value of 4.4 mBqm -3 , we estimate that this air was 
diluted by a factor of 10 5 -10 6 prior to reaching the bound- 
ary of the continental U.S. This level of dilution is not 
surprising given the transport patterns mentioned above. 
From the presence and/or absence of certain fission iso- 



topes we can draw several conclusions on their origin: (i) 
The value of the ratio of 134 Cs to 137 Cs activities of ~0.7 is 
indicative of the release of the fission products from a nu- 



clear reactor and not from nuclear weapons (Devell et al. 



1986). (ii) The presence of the relatively short lived iso- 



topes 131 I and 132 Te shows that the fission products had 
been released primarily from recently active fuel rods as 
opposed to spent fuel, (iii) The notable absence of 20.8-h 
133 I in our spectra, together with the known steady state 
activity ratio of 133 I to 131 I of -2 (|Devell et a~ 



1986 



Iso- 



tope Project 2011) allows us to put a lower limit on the 



time between the end of steady state nuclear fuel burning 
and the arrival of the fission products at our location. A 
lower limit of 5.7 days (95% C.L.) was derived from the 
PH2 measurement to be compared with 6.6 days delay 
between the earthquake and the end of PH2 filter expo- 
sure. The PH4 measurement resulted in a lower limit of 
8.1 days (95% C.L.) for the end of exposure of 8.4 days af- 
ter the earthquake. Given the modelled transport time of 
5-6 days from above this means that the reactor success- 
fully shut down at the time of the earthquake and that 
nuclear reactions were largely stopped at the time of the 
release of the fission products. Also we found a small indi- 
cation for a 133 I peak in our PH2 data set. The peak with 
22±8 counts/hour per day of exposure was not significant 
enough to pass our 3tr cut. The reactor shut down time 
derived from the insignificant 133 I peak is 6.8±0.5 days 
before the end of PH2 exposure which is statistically com- 
patible with the time of the earthquake, (iv) It is striking 
that we see only three of the many possible fission product 
elements. This points to a specific process of release into 
the atmosphere. The exact process and why it would be se- 
lective requires further investigation, but we can speculate 
that the release of fission products to the atmosphere is the 
result of evaporation of contaminated steam, in which, e.g., 
Csl is very soluble. Chernobyl debris, conversely, showed 
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NOAAHYSPLIT MODEL 
Forward trajectories starting at 1000 UTC 12 Mar 1 1 
GDAS Meteorological Data 




This is not a NOAA product. It was produced by a web user. 

Job ID: 307297 Job Start: Thu Apr 717:17:21 UTC2011 

Source 1 [at.: 37.75 Ion.: 140.43 hgts: 500, 1 000, 1 500 m AGL 

Trajectory Direction: Forward Duration: 240 hrs 

Vertical Motion Calculation Method: Model Vertical Velocity 

Meleoroloay: OOOOZoa Mar 2011 -GDAS1 



(a) 




(b) 



Figure 5: (Colour online) Calculated trajectories for airmasses re- 
leased at the site of the Fukushima Dai-ichi reactor (a) and the 
850 mb geopotential heights, wind vectors and temperatures for 17 
March 2011, 12 UTC (b). For more details see text. 



a much broader spectrum of elements ( Devell et al. 1986 ) 



reflecting the direct dispersal of active fuel elements. 



4. Conclusions 

We measured the arrival of airborne fission products 
from the Fukushima Dai-ichi, Japan nuclear reactor inci- 
dent in Seattle, WA, USA. The first fission products ar- 
rived between 17-18 March about 7 days after the earth- 
quake and in agreement with other reported detections 



of radionuclides in the western United States (RadNet 



2011 



2011 



|Taylor[ [20TT| |Norman etH] [20TT] |Bowyer et ah 



Our models of the transport of air masses across 
the Pacific point to a typical transport time of 5 to 6 days. 
This agrees with the first reported explosion at Fukushima 



one day after the earthquake (Tokyo Electric Power Com- 



pany TEPCO 20111, augmented by additional contribu- 



tions from subsequent explosions. 

The detected fission isotopes are 131 1, 132 I, 132 Te, 134 Cs, 
and 137 Cs with the highest activity observed for 131 I of 
4.4±1.3 mBqm~ 3 on 19-20 March. Our measurements 
were only sensitive to radionuclides attached to particu- 
late mater and were thus not sensitive to the part of 131 I 
that is transported in its gaseous form. The presence of 
the aforementioned isotopes clearly points to the release of 
the fission products from recently active fuel elements. At 
the same time, the absence of the short lived isotope 133 I 
indicates that nuclear reactions must have been stopped 
successfully at the time of the earthquake as we would 
otherwise have detected its activity. 

This set of measurements over a period of 23 days 
provide a quantitative basis for further modeling on how 
radioactive fission products are transported in the atmo- 
sphere and will hopefully be included in a global analysis. 
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